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Abstract. We prove quaternionic contact versions of two of Obata's sphere theorems. On a compact 
quaternionic contact (qc) manifold of dimension bigger than seven and satisfying a Lichnerowicz type lower 
bound estimate we show that if the first positive eigenvalue of the sub-Laplacian takes the smallest possible 
value then, up to a homothety of the qc structure, the manifold is qc equivalent to the standard 3-Sasakian 
sphere. The same conclusion is shown to hold on a non-compact qc manifold which is complete with 
respect to the associated Riemannian metric assuming the existence of a function with traceless horizontal 
Hessian. The third result of the paper is a qc version of Liouville's theorem showing that a qc-conformal 
diffcomorphism between open connected sets of the 3-Sasakian sphere is a restriction of an element of the 
qc-conformal automorphism group of the sphere. 
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1. Introduction 

Motivated by the classical Lichnerowicz [62] and Obata [71] theorems, earlier papers of the authors 
[40, 41] established a Lichnerowcz type lower bound estimate for the first eigenvalue of the sub-Laplacian 
on a compact quatcrnionic contact (qc) manifold. The case of equality in the lower bound estimate (Obata- 
type theorem) was settled in the special case of a 3-Sasakian compact manifold where it was shown that the 
lower bound for the first eigenvalue of the sub-Laplacian is achieved if and only if the 3-Sasakian manifold 
is isometric to the standard 3-Sasakian sphere. Quatcrnionic contact (qc) structures were introduced by O. 
Biquard [6] and arc modeled on the conformal boundary at infinity of the quaternionic hyperbolic space. 
Thus, manifolds equipped with a qc structure are examples of sub-Riemannian geometries. The (locally) 
3-Sasakian manifolds were characterized in [37, 43] by the vanishing of the torsion tensor of the Biquard 
connection. The qc geometry was a crucial geometric tool in finding the extremals and the best constant in 
the L 2 Folland-Stcin Sobolcv-type embedding, [27, 28], completely described on the quaternionic Heisenberg 
groups, [38, 39]. 

In this paper we prove the full qc version of Obata's results for a general qc manifold of dimension bigger 
than seven. We find that the equality case of Lichnerowcz' type inequality on a compact qc manifold of 
dimension at least eleven can be achieved only on the 3-Sasakian spheres. More general, we show that on 
a complete with respect to the associated Ricmannian metric qc manifold a certain (horizontal) Hessian 
equation, cf. (1.6), allows a non-trivial solution if and only if the manifold is qc nomothetic to the standard 
3-Sasakian sphere. 

The qc seven dimensional case was considered in [41], however, the general qc Obata results in dimension 
seven remain open. 

Turning to some details, let us recall the mentioned classical results. Using the classical Bochner- 
Weitzenbock formula Lichnerowcz [62] showed that on a compact Ricmannian manifold (M, h) of dimension 
n for which the Ricci curvature is greater than or equal to that of the round unit n-dimensional sphere S n (l), 
i.e., Ric(X, Y) > (n — l)h(X, Y) the first positive eigenvalue Ai of the (positive) Laplace operator is greater 
than or equal to the first eigenvalue of the sphere, Ai > n. Subsequently Obata [71] proved that equality is 
achieved if and only if the Ricmannian manifold is isometric to S n (l) by noting that the trace- free part of 
the Riemannian Hessian of an eigenfunction / with eigenvalue A = n vanishes, i.e., it satisfies the system 

(1.1) {V h ?f = -fh 

after which he defines an isomctry using analysis based on the geodesies and Hessian comparison of the 
distance function from a point. 

In fact, Obata showed that on a complete Ricmannian manifold (M,h) equation (1.1) allows a non- 
constant solution if and only if the manifold is isometric to the unit sphere. In this case, the eigenfunctions 
corresponding to the first eigenvalue are the solutions of (1.1). Later, Gallot [30] generalized these results 
to statements involving the higher eigenvalues and corresponding eigenfunctions of the Laplace operator. 

The interest in relations between the spectrum of the Laplacian and geometric quantities justified the 
interest in Lichnerowicz-Obata type theorems in other geometric settings such as Riemannian foliations (and 
the eigenvalues of the basic Laplacian) [58, 57], [50] and [72], to CR geometry (and the eigenvalues of the 
sub-Laplacian) [34], [4], [17, 15, 16], [18], [20], [59], and to general sub-Riemannian geometries, see [5] and 
[36]. In the CR case, Greenleaf [34] gave a version of Lichnerowicz' result showing that if a compact strongly 
pseudo-convex CR manifold M of dimension 2n + 1, n > 3 satisfies a Lichnerowicz type inequality 

Ric(X, Y) + 4A(X, JY) >(n + l)g(X, X) 

for all horizontal vectors X, where Ric and A are, correspondingly, the Ricci curvature and the Webster 
torsion of the Tanaka- Webster connection (in the notation from [47, 44]), then the first positive eigenvalue 
Ai of the sub-Laplacian satisfies the inequality Ai > n. The standard (Sasakian) CR structure on the 
sphere achieves equality in this inequality. Following [34] the above cited results on a compact CR manifold 
focused on adding a corresponding inequality for n = 1, 2 or characterizing the equality case mainly in the 
vanishing Webster-torsion case (the Sasakian case). The general case on a compact CR manifold satisfying 
the Lichnerowicz type condition was proved in [60, 61] using the results and the method of [45] while 
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introducing a new integration by parts step proving the vanishing of the Webster torsion assuming the first 
eigenvalue is equal to n (for the three dimensional case see [46]). On the other hand, a generalization of the 
Obata result in the complete non-compact case was achieved in [45] , where the standard Sasakian structure 
on the unit sphere was characterized through the existence of a non-trivial solution of a (horizontal) Hessian 
equation on a complete with respect to the associated Riemannian metric CR manifold with a divergence 
free Webster torsion. To the best of our knowledge the case of a general torsion remains still open. 

The main purpose of this paper is to prove the qc version of both results of Obata under no extra 
assumptions on the Biquard' torsion when the dimension of the qc manifold is at least eleven, cf. Theorem 1.2 
and Theorem 1.3. 

The quatcrnionic contact version of the Lichnerowicz' result was found in [40] in dimensions grater than 
seven and in [41] in the seven dimensional case. The following result of [40] gives a lower bound on the 
positive eigenvalues of the sub-Laplacian on a qc manifold. 

Theorem 1.1 ([40]). Let (M,rj, g,Q) be a compact quaternionic contact manifold of dimension 4n + 3 > 7. 
Suppose that there is a positive constant kg such that the qc Ricci tensor and torsion of the Biquard connection 
satisfy the inequality 

, N , , 2(4n + 5) n/ , 6(2n 2 + 5n-l) , , 

(1.2) Ric(X,X) + -L—Jt^X) + A___2 W > kog{X ,X). 

Then any eigenvalue X of the sub-Laplacian A satisfies the inequality 

n 

A> — 

n + 2 

The equality case of Theorem 1.1 is achieved on the 3-Sasakian sphere. It was shown in [39], see also [3], 
that the eigenspace of the first non-zero eigenvalue of the sub-Laplacian on the unit 3-Sasakian sphere in 
Euclidean space is given by the restrictions to the sphere of all linear functions. 

The main results of this paper are the following three theorems. 

Theorem 1.2. Let (AT, 77,(7, Q) be a compact quaternionic contact manifold of dimension 4n + 3 > 7 whose 
qc-Ricci tensor and torsion of the Biquard connection satisfy the inequality (1.2). Then, the first positive 
eigenvalue X of the sub-Laplacian A satisfies the equality 

(1.3) X = -^-k 

n + 2 

if and only if the qc manifold (M, g, Q) is qc-homothetic to the unit Hn+ 3) -dimensional 3-Sasakian sphere. 

According to [40, Remark 4.1], under the conditions of Theorem 1.1, an eigenfunction / corresponding 
to the the first non-zero eigenvalue as in (1.3), A/ = ^p^ko f satisfies a linear PDE system, namely, the 
horizontal Hessian of / is given by (see Corollary 4.2 in the Appendix) 

1 3 
(1-4) Vdf(X,Y) = - - k fg(X, Y)-J2 Y). 

^ ' s— 1 

This brings us to our second main result. 

Theorem 1.3. Let (M, 77, g, Q) be a quaternionic contact manifold of dimension 4n+3 > 7 which is complete 
with respect to the associated Riemannian metric 

(1-5) h = g+(i ll f + {i l2 f + { m f. 

Suppose there exists a non-constant smooth function f whose horizontal Hessian satisfies 

3 

(1-6) Vd/(X, Y) = -fg(X, Y)-J2 df(£ s )u s (X, Y). 

Then the qc manifold (M, rj,g, Q) is qc homothetic to the unit (4n+ 3) -dimensional 3-Sasakian sphere. 
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Clearly Theorem 1.3 implies Theorem 1.2 since any Riemannian metric on a compact manifold is complete 
and a qc-homothcty, cf. Definition 3.12, allows us to reduce to the case fco = 4(n + 2), which turns (1.4) in 
(1.6). 

We achieve the proof of Theorem 1.3 by showing first that M is isometric to the unit sphere S in+3 and 
then that M is qc-cquivalent to the standard 3-Sasakian structure on S 4n+3 . To this effect we show that the 
torsion of the Biquard connection vanishes and in this case the Riemannian Hessian satisfies (1.1) after which 
we invoke the classical Obata theorem showing that M is isometric to the unit sphere. In order to prove the 
qc-equivalence part we show that the qc-conformal curvature vanishes, which gives the local qc conformal 
equivalence with the 3-Sasakian sphere due to [42, Theorem 1.3], and then use the below Liouville-type 
result Theorem 1.4, which implies the existence of a global qc-conformal map between M and the 3-Sasakian 
sphere, cf. sub-section 3.8. 

Theorem 1.4. Let £ C S 4n+3 be a connected open subset of S 4n+3 . If F : £ — > S 4n+3 is a qc-conformal 
transformation then F is the restriction to £ of an element of PSp(n +1,1) - the isometry group of the 
quaternionic hyperbolic space. 

The proof of Theorem 1.4 and background to this result can be found in sub-sections 3.8 and 3.9. 

A version of Theorem 1.1 when n = 1 was established in [41, Theorem 1.1] assuming the positivity of the 
P-function of any eigcnfunction. In the Appendix, for completeness, we recall the notion of the P-function 
introduced in [41] and give a different proof of Theorem 1.1 based on the positivity of the P-function in the 
case n > 1 established in [41, Theorem 3.3]. As a corollary of the proof, we show the validity of (1.4) for 
any eigcnfunction of the sub-Laplacian with eigenvalue given by (1.3). 

Convention 1.5. 

a) We shall use X, Y, Z, U to denote horizontal vector fields, i.e. X,Y, Z,U £ H . 

b) {ei, . . . , e4„} denotes a local orthonormal basis of the horizontal space H. 

c) The summation convention over repeated vectors from the basis {ei, . . . , e^n} will be used. For example, 
for a (0,4)-tensor P, the formula k = P(e&, e a , e a , e&) means k = X^ ™ b=1 P( e b, e a, e a , e&)- 

d) The triple (i,j,k) denotes any cyclic permutation of (1,2,3). 

e) The sum J^fyfc) means the cyclic sum. For example, 

dfihX^Y, Z) = df{hX)uj 2 (Y, Z) + df(I 2 X)uj 3 (Y, Z) + df{I^X)w\ (Y, Z). 

(ijk) 

e) s will be any number from the set {1,2,3}, s £ {1,2,3}. 

Acknowledgments The research is partially supported by the Contract "Idci", DID 02-39/21.12.2009. 
S.I and A. P. are partially supported by the Contract 130/2012 with the University of Sofia 'St.Kl.Ohridski'. 
D.V. would like to thank Professor Luca Capogna for some useful comments. 

2. Quaternionic contact manifolds 

In this section we will briefly review the basic notions of quaternionic contact geometry and recall some 
results from [6], [37] and [42] which we will use in this paper. 

It is well known that the sphere at infinity of a non-compact symmetric space M of rank one carries a 
natural Carnot-Caratheodory structure, see [66, 70]. In the real hyperbolic case one obtains the conformal 
class of the round metric on the sphere. In the remaining cases, each of the complex, quaternion and 
octonionic hyperbolic metrics on the unit ball induces a Carnot-Caratheodory structure on the unit sphere. 
This defines a conformal structure on a sub-bundle of the tangent bundle of co-dimension dimK — 1, where 
K = C, H, O. In the complex case the obtained geometry is the well studied standard CR structure on the 
unit sphere in complex space. Quaternionic contact (qc) structure were introduced by O. Biquard, sec [6], and 
are modeled on the conformal boundary at infinity of the quaternionic hyperbolic space. Biquard showed 
that the infinite dimensional family [55] of complete quaternionic-Kahler deformations of the quaternion 
hyperbolic metric have conformal infinities which provide an infinite dimensional family of examples of qc 
structures. Conversely, according to [6] every real analytic qc structure on a manifold M of dimension at 
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least eleven is the conformal infinity of a unique quaternionic-Kahler metric denned in a neighborhood of M. 
Furthermore, [6] considered CR and qc structures as boundaries of infinity of Einstein metrics rather than 
only as boundaries at infinity of Kahler-Einstein and quaternionic-Kahler metrics, respectively. In fact, in 
[6] it was shown that in each of the three cases (complex, quaternionic, octoninoic) any small perturbation 
of the standard Carnot-Carathcodory structure on the boundary is the conformal infinity of an essentially 
unique Einstein metric on the unit ball, which is asymptotically symmetric. In the Riemannian case the 
corresponding question was posed in [26] and the perturbation result was proven in [33] . 

Another natural extension of an interesting Riemannian problem is the quaternionic contact Yamabe 
problem, a particular case of which [31, 79, 37, 38] amounts to finding the best constant in the L 2 Folland- 
Stcin Sobolev-type embedding and the functions for which the equality is achieved, [27] and [28] with a 
complete solution on the quaternionic Heisenberg groups given in [38, 39]. 

2.1. Quaternionic contact structures and the Biquard connection. A quaternionic contact (qc) man- 
ifold (M, rj, g, Q) is a An + 3-dimcnsional manifold M with a codimension three distribution H locally given 
as the kernel of a 1-form r\ = (771, 772, 773) with values in K 3 . In addition H has an Sp(n)Sp(l) structure, that 
is, it is equipped with a Riemannian metric g and a rank-three bundle Q consisting of endomorphisms of H 
locally generated by three almost complex structures I\ , I2 , ^3 on H satisfying the identities of the imagi- 
nary unit quaternions, I\Ii = = I3, I1I2I3 = — which are hermitian compatible with the metric 
g(I s -, I s -) = <?(■, ■) and the following compatibility condition holds 2g(I s X, Y) — di] s (X, Y), X,Y € H. 

The transformations preserving a given quaternionic contact structure 77, i.e., fj = u^n for a positive 
smooth function u and an 5*0(3) matrix '5 with smooth functions as entries are called quaternionic contact 
conformal (qc- conformal) transformations, see Deinition 3.12. If the function /1 is constant 77 is called qc- 
homothetic to rj. The qc conformal curvature tensor W qc , introduced in [42], is the obstruction for a qc 
structure to be locally qc conformal to the standard 3-Sasakian structure on the (4n + 3)-dimcnsional sphere 



A special phenomena, noted in [6], is that the contact form rj determines the quaternionic structure and 
the metric on the horizontal distribution in a unique way. 

On a qc manifold with a fixed metric g on H there exists a canonical connection defined first by O. 
Biquard in [6] when the dimension (4n + 3) > 7, and in [25] for the 7-dimensional case. Biquard showed that 
there is a unique connection V with torsion T and a unique supplementary subspace V to H in TM, such 
that: 

(i) V preserves the decomposition H © V and the Sp(n)Sp(l) structure on H, i.e. Vg = 0, Vcr £ T(Q) for 
a section a £ T(Q), and its torsion on H is given by T(X, Y) = —[X, Y\\v) 

(ii) for £ £ V, the endomorphism T(£, .)i H of H lies in (sp(n) © sp(l)) ± C gl(An); 

(iii) the connection on V is induced by the natural identification (p of V with the subspace sp(l) of the 
endomorphisms of H, i.e. V</? = 0. 

This canonical connection is also known as the Biquard connection.MVhen the dimension of M is at least 
eleven [6] also described the supplementary distribution V, which is (locally) generated by the so called 
Reeb vector fields {^1,^2,^3} determined by 



where j denotes the interior multiplication. If the dimension of M is seven Duchemin shows in [25] that if we 
assume, in addition, the existence of Reeb vector fields as in (2.1), then the Biquard result holds. Henceforth, 
by a qc structure in dimension 7 we shall mean a qc structure satisfying (2.1). 

Notice that equations (2.1) are invariant under the natural 5(9(3) action. Using the triple of Reeb vector 
fields we extend the metric g on H to a metric h on M by requiring span{^\, £2, £3} = V _L H and h(^ s ,^k) = 
S s k- The Riemannian metric h as well as the Biquard connection do not depend on the action of 50(3) on 
V, but both change if 77 is multiplied by a conformal factor [37]. Clearly, the Biquard connection preserves 
the Riemannian metric on TM, V/i = 0. Since the Biquard connection is metric it is connected with the 



[37, 42]. 
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>7s(fk) = $sk, 



(£ s _id77 s )|if = 0, 
= -(6fcjd?7 s )|H, 
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Levi-Civita connection V h of the metric h by the general formula 

(2.2) h(V A B, C) = h(W h A B, C) + - h(T(A, B),C)- h(T(B, C), A) + h(T(C, A),B) , A,B,C 6 T(TM). 

The covariant derivative of the qc structure with respect to the Biquard connection and the covariant 
derivative of the distribution V are given by 

(2.3) VJi = —<Xj ® 4 +aik ® Ij, V& = -ctj ® & + a-k ® £?• 

The vanishing of the sp(l)-connection 1-forms on H implies the vanishing of the torsion endomorphism of 
the Biquard connection (see [37]). 

The fundamental 2-forms to s of the quaternionic structure Q arc defined by 

(2.4) 2lu s]h = dn s \ H , £jw s = 0, £eV. 
Due to (2.4), the torsion restricted to H has the form 

(2.5) T(X, Y) = -[X, Y] lv = 2wi(X, Y)& + 2u 2 (X, Y)& + 2uj 3 (X, Y)&. 

2.2. Invariant decompositions. An endomorphism \& of H can be decomposed with respect to the quater- 
nionic structure (Q, g) uniquely into four 5p(n)-invariant parts ^> = \E»+++ + v|> H + <J/ 1 + <J/ ^, where 
vf+++ commutes with all three h, commutes with I\ and anti-commutes with the others two and etc. 

The two 5p(n)5p(l)-invariant components ^[ 3 j = ^[-i] = ^ + ^ 1 + ^ h ar e determined 

by 

* = *[ 3 ] <^=> 3* + il^/l + / 2 *-^2 + h^h = 0, 

* = * - Ii^/i - h^h - h^h = 0. 

With a short calculation one sees that the S'p(n)S'p(l)-invariant components are the projections on the 
cigenspaces of the Casimir operator T = I-y ® I\ + I 2 ® I2 + -^3 ® I3 corresponding, respectively, to 
the eigenvalues 3 and —1, see [12]. If n = 1 then the space of symmetric endomorphisms commuting with 
all I s is 1-dimcnsional, i.e. the [3]-component of any symmetric endomorphism A> on H is proportional to 
the identity, ^3 = — ^-Idm- Note here that each of the three 2-forms w s belongs to its [-l]-component, 
u s = w a r_i] and constitute a basis of the Lie algebra sp(l). 



2.3. The torsion tensor. The properties of the Biquard connection are encoded in the properties of the 
torsion endomorphism 7^ = T(£, •) : H — >• H, £ 6 V. Decomposing the endomorphism S (sp(n) + 
sp(l))' L into its symmetric part T® and skew-symmetric part b^,T^ = T® + b$, O. Biquard shows in [6] that 
the torsion XJr is completely trace- free, trT^ = trT^ o I a = 0, its symmetric part has the properties T9.Ii = 
-UTl h{TD+- = /i(T°)-+-, h{Tl)~+- = / 2 (T°)--+, ^ h{Tl)-+ = h(Tl )+~, where the 

upperscript H — | — h means commuting with all three 7,, H — indicates commuting with ii and anti- 

commuting with the other two and etc. The skew-symmetric part can be represented as b^ t = IiU, where u 
is a traceless symmetric (1,1) -tensor on H which commutes with Ii,l2,H. Therefore we haveT^ = T®. + IiU. 
If n = 1 then the tensor u vanishes identically, u = 0, and the torsion is a symmetric tensor, = T®. 

Any 3-Sasakian manifold has zero torsion endomorphism, = 0, and the converse is true if in addition 
the qc scalar curvature (see (2.6)) is a positive constant [37] (the case of negative qc-scalar curvature can 
be treated very similarly, see [43, 44]. We remind that a (An + 3)-dimcnsional Ricmannian manifold (M, g) 
is called 3-Sasakian if the cone metric g c = t 2 h + dt 2 on C = M x R + is a hyper Kahler metric, namely, 
it has holonomy contained in Spin + 1) [9]. A 3-Sasakian manifold of dimension (An + 3) is Einstein with 
positive Ricmannian scalar curvature (An + 2)(4n + 3) [51] and if complete it is a compact manifold with a 
finite fundamental group (see [8] for a nice overview of 3-Sasakian spaces). 

2.4. Torsion and curvature. Let R = [V, V] — V[ ] be the curvature tensor of V and the dimension is 
An + 3. We denote the curvature tensor of type (0,4) and the torsion tensor of type (0,3) by the same letter, 
R(A, B, C, D) := h(R(A, B)C, D), T(A 7 B, C) := h(T(A, B),C), A 7 B 7 C,D e T(TM). The Ricci tensor, 
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the normalized scalar curvature, the Ricci 2-forms and the Ricci type-tensor £ of the Biquard connection, 
called qc-Ricci tensor Ric, normalized qc-scalar curvature 5, qc- Ricci forms p s , respectively, are given by 



Ric(A, B) = R(e b , A, B, e b ), 8n(n + 2)5 = R(e b , e a , e a ,e b ), 



(2.6) 



Ps(A, B) 



1 



R(A,B,e a ,I s e a ), C S (A,B) 



1 



R(e a ,A,B,I s e a ). 



An in 
The sp(l)-part of R is determined by the Ricci 2-forms and the connection 1-forms by 

(2.7) R(A, B,£i,£j) = 2p k (A, B) = (da k + a, t A atj)(A,B), A,B £ T(TM). 

The two 5p(n)5p(l)-invariant tracc-frec symmetric 2-tcnsors T°(X,Y) = g^T^h + T°J 2 +Tgl s )X,Y), 
U(X,Y) = g(uX,Y) on H, introduced in [37], have the properties: 



(2.8) 



T°{X, Y) + T°(hX, hY) + T°(I 2 X, I 2 Y) + T°{I 3 X, I 3 Y) = 0, 
U(X,Y) = U(hX,hY) = U(hX,I 2 Y) = U{hX,hY). 



In dimension seven (n — 1), the tensor U vanishes identically, [7 = 0. 

We shall need the following identity taken from [42, Proposition 2.3] 4T°(^ S , I S X, Y) 
T°(I S X, I S Y) which implies the formula 

(2.9) T(&, I S X, Y) = T°(£ S ,I S X, Y) + g(I s uI s X, Y) = \ [t°(X, Y) - T°{I S X, I S Y) 



T°(X,Y) 



U(X,Y). 



We recall that a qc structure is said to be qc-Einstein if the horizontal qc-Ricci tensor is a scalar multiple of 
the metric, Ric(X, Y) = 2(n + 2)Sg(X,Y). The horizontal Ricci-type tensor can be expressed in terms of 
the torsion of the Biquard connection [37] (see also [38, 42]). We collect below the necessary facts from [37, 
Theorem 1.3, Theorem 3.12, Corollary 3.14, Proposition 4.3 and Proposition 4.4] with slight modification 
presented in [42] 

Ric{X, Y) = (2n + 2)T°(X, Y) + {An + IQ)U{X, Y) + 2(n + 2)Sg(X, Y), 
p s (X, I S Y) = -\[t°(X, Y) + T a (I s X, I,Y)] - 2U(X, Y) - Sg(X, Y), 



(2.10) 



UX,I S Y) 



2n + 1 n , % In, 

— r T\X,Y) + -T\l s XJ s Y) 



2 -^U(X,Y )+ S - 9{ X, n 



T(&,fc) = -ss; k - s = -/i(T(a,6),6), 

g(T(^),X) = -pkihX^) = -p^X,^) = -h([^,^},X). 

For n = 1 the above formulas hold with U = 0. Hence, the qc-Einstein condition is equivalent to the 
vanishing of the torsion endomorphism of the Biquard connection. In this case the normalized qc scalar 
curvature S is constant and the vertical distribution V is integrable provided n > 1. If 5 > then the qc 
manifold is locally 3-Sasakian [37], (sec [43] for the negative qc scalar curvature). 
We shall also need the general formula for the curvature [42, 44] 



(2.11) R(&,X,Y,Z) 
1 -, 
~ 4 



■(y x U)(IiY, Z) + Uj (X, Y) Pk {I i Z, e 4 ) - u k (X, Y)p j {I i Z, £<) 
(\7 Y T°)(hZ,X) + (V Y T ){ZJ t X) \ + -\(V z T°){hY,X) + (V Z T°)(Y, hX) 



- Uj (X, Z) Pk (I l Y, &) + u k (X, Z)p j {I i Y, fc) - 0Jj(Y, Z)p k (I i X, &) + Wfc (Y, Z)p j {I i X, 
where the Ricci two forms are given by, cf. [42, Theorem 3.1] or [44, Theorem4.3.11] 

6(2n + l)p a (£.,X) = (2n + l)X(S) + i(V eo T°)[(e ,X) - 3(I s e a , I.X)\ - 2(V ea U)(e a ,X), 

( 2 ' 12) 6(2n + l)pi{ij, I k X) = (2n - l)(2n + l)X(S) - i(V ea T°)[(4n + l)(e a , X) + 3(/ 4 e a , I t X)] 

-A(n + l)(S7 ea U)(e a ,X). 
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2.5. The Ricci identities, the divergence theorem. We shall use repeatedly the following Ricci iden- 
tities of order two and three, see also [42] and[40]. Let £j, i = 1, 2,3 be the Reeb vector fields, / a smooth 
function on the qc manifold M and V/ its horizontal gradient, g(Vf, X) = df(X). We have: 

3 

v 2 /(A,y) - v 2 /(y,x) = -2^> s (x,y)d/(£ s ), 

s=l 

V 2 /(A, 6) - V 2 /(6, X) = T(£ a ,X, Vf), 

3 

R(X, Y, Z, Vf) - 2 Wa(X, F)V 2 /(6, 

s=l 

3 

2rf/(O fe (x, y) + 2df^ k ) Pj (x, y) - 2 £ Ws (x, r)v 2 (6, 6), 

8=1 

-/*(&, x, y, v/) - v 2 /(T(e s , x), y), 

-V 2 / (T (6, X) , y) - V 2 / (X, T (£„ Y)) - df ((V X T) 

-R(£ s ,X,Y,Vf). 

The horizontal sub-Laplacian A/ and the norm of the horizontal gradient V/ of a smooth function / on 
M are defined respectively by 

Af = -tr a H {V 2 f) = V*df = - V 2 /(e a ,e a ), |V/| 2 = df (e a ) df (e a ) . 

The function / is an cigcnfunction with eigenvalue A of the sub-Laplacian if, for some constant A we have 

(2.14) A/ = A/. 

From the Ricci identities we have the following formulas for the traces through the almost complex structures 
of the Hessian 

(2.15) g(V 2 f, «,) = V 2 /(e a , 7 s e Q ) = -4ru$f (&)■ 

For a fixed local 1-form 77 and a fix s £ {1, 2, 3} the form Vol v = t}\ A 772 A 773 Aw 2 ™ is a locally defined volume 
form. Note that Vol v is independent on s as well as it is independent on the local one forms 771, 772, 773. Hence 
it is globally defined volume form denoted with Vol v . The (horizontal) divergence of a horizontal vector 
field/one- form cr € A 1 (if) defined by V* a = —tr\iiVcr = — Vcr(e a , e a ) supplies the "integration by parts" 
formula [37], see also [79], 

(2.16) / (VV) Vo^ = 0. 

Jm 

3. Proof of the main Theorems 

The proof of Theorem 1.3 is lengthy and requires a number of steps which we present in the following 
sub-sections. Throughout this section we shall assume the assumptions of Theorem 1.3. In particular, / is a 
non-constant smooth function whose horizontal Hessian satisfies (1.6). Our first step is to show the vanishing 
of the torsion tensor, T° = and U = 0. We start by expressing the remaining parts of the Hessian (w.r.t. 
the Biquard connection) in terms of the torsion tensors and show that / satisfies an elliptic equation on M. 
A simple argument shows that T°(7 S V/, V/) = U(I s Vf,Vf) = 0, s = 1,2,3. Furthermore, using the [— 1]- 
component of the curvature tensor we show that T°(I S Vf, ItVf) = 0, s, t £ {1, 2, 3}, s ^ t. In addition, we 
determine the torsion tensors T° and U in terms of the horizontal gradient of / and the tensor U(Vf, Vf). 
The analysis proceeds by finding formulas of the same type for the covariant derivatives of T° and U. Thus, 
the crux of the matter in showing that the torsion vanishes is the proof that U(Vf, Vf) = 0. This fact 
will be achieved with the help of the Ricci identities, the contracted Bianchi second identity and thus far 
established results. In the next step of the proof of Theorem 1.3 we compute the Riemannian Hessian of /, 
with respect to the Levi-Civita connection of the metric (1.5) which allow us to invoke Obata's result thus 
proving that M equipped with the Riemannian metric (1.5) is homothetic to the unit sphere in quaternion 
space. The final step is to show that M is qc-homothetic to the (4n + 3)-dimcnsional 3-Sasakian unit sphere. 



V 3 f(X,Y,Z) 

(2.13) 

V 3 f(X,Y,^) 

v 3 /(6,A,y) 

V 3 f(t s ,X,Y) 



V 3 f(Y,X,Z) = - 

V 3 f(Y,X,^) = - 

V 3 f(X,i„Y) = - 
V 3 f(X,Y,£ s ) = - 
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Here, we employ a standard monodromy argument showing that a compact simply connected locally qc- 
conformally flat manifold is globally qc-conformal to the 3-Sasakian unit sphere. The key is Theorem 1.4, 
which is a generalization of the Liouville theorem, showing that every qc-conformal transformation between 
open subsets of the 3-Sasakian unit sphere is the restriction of a global qc-conformal transformation, i.e., an 
element of the group PSp(n + 1, 1), sec subsection 3.8 for further details. 

3.1. Some basic identities. We start our analysis by finding a formula for the the third covariant derivative 
of a function which satisfies (1.6). 

Lemma 3.1. With the assumptions of Theorem 1.3 we have the following formula for the third covariant 
derivative of the function f , 

3 

(3.1) V 3 f(A,X,Y) = -df(A)g(X,Y)-Y,Us(X,Y)V 2 f(A,Z s ), A e T(TM). 

s=l 

Proof. The claimed formula is obtained by differentiating the Hessian equation (1.6). Indeed, the covariant 
derivative along A S T(TM) of (1.6) gives 



V 3 f(A,X,Y) = ~df(A)g(X,Y) 



[V 2 /(A Y) + df(V A Z s )oj s {X, Y) + df{£ s ) (\7 a cj s ) (X, Y)] 



s=l 



which together with (2.3) gives the identity, cf. also Convention 1.4 e), 

V 3 f{A 7 X, Y) = -df(A)g(X, Y) - £ [V 2 f{A, Y) + #(V A &)<*(*. Y ) + 4f (6) ( v ^) (X, Y)] 

(ijk) 

3 

= -df(A)g(X,Y) -J2 [V 2 f(A,Z t )uj t (X,Y)] 
t=i 

- [-ajW&l + ^W^^i^y)- [-^(A)MXX) + MAH(x,Y)]df(^) 

(ijk) (ijk) 

3 

= -rf/(%(i,F)-^[v 2 /(A6K(x,y)] , 
*=i 

which completes the proof. □ 

After this technical Lemma, our first goal is to find a formula for the curvature tensor R(Z,X,Y,Vf), 
for / satisfying (1.6), using Lemma 3.1 with A = Z, the Rlcci identities (2.13), and the properties of the 
torsion. In fact, after some standard calculations it follows 

(3.2) R(Z, X, Y, Vf) = [df{Z)g{X, Y) - df{X)g{Z, Y) 

3 

+ J2 [vdf(Z s ,Z)u s (X,Y) -Vdf(Z s ,X)u s (Z,Y) -2Vdf(£ s ,Y)u s (Z,X) 

s=l 

3 

+ Y, [t(Zs,Z, Vf)uj s (X, Y) - T(&, X, Vf)uj s (Z, Y) . 

s=l 

By taking traces in (3.2) we can derive formulas for the various contracted tensors (2.6). We shall use the 
following, 

3 3 

Ric(Z, V/) = (4n -l)df(Z)-Y T{i s ,I s Z, V/) - 3^ Vdf(£ s ,I s Z), 

(3.3) 8=1 s=1 

4n&(IiZ, V/) = -df(Z) + {An - 1)T(&, hZ, V/) + TfaljZ, Vf) + T(&, I k Z, V/) 

+ {An + l)Wdf{^hZ) - VdffaJjZ) - V4f(&, hZ). 
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The above formulas imply some other basic identities to which we turn next. Note that with the help of 
(2.10) we can rewrite the Lichnerowicz type assumption (1.2) in the form 

(3.4) L(X, X) d = f 2(ti + 2)Sg(X, X) + a' n T°(X, X) + (3' n (X, X) > k Q g(X, X), X £ H, 

, _ 2(2n + 3)(n + 2) _ 4(2n - l)(n + 2) 2 

2n + l ' Pn ~ (2n+l)(n-l) ' 

which allows to write the first claim of the following Lemma in the form L(Z, V/) = for all Z £ H whenever 
/ satisfies (1.6) taking fco = 4(n + 2). 

Lemma 3.2. With the assumptions of Theorem 1.3, the next identity holds true 

(3.5, (5 - we, + £±5n* v/) + 2 %:+Z-t u{z ' v/) - °- 

Furthermore, we have 

(3.6) T°(I s Wf, V/) - 0, U(I s Vf, V/) = 0. 
Proof. The first equations in (3.3) and (2.10) together with (2.9) imply 

3 

(3.7) 3 Vdf(£ s ,I s Z) = An - 1 - (2n + 4)5 df{Z) - (2n + 3)T°(Z, V/) - (An + 7)U(Z, V/). 

s=l 

The sum over 1, 2, 3 of the second equality in (3.3) together with the third equality of (2.10) and (2.9) gives 

3 

(3.8) (An - 1) ^ Vdf(£ s ,I s Z) = (3 - 6nS)df(Z) - (2n + 3)T°)(Z, V/) - 3(7(Z, V/). 

s=l 

Subtracting (3.7) from (3.8) we obtain 

3 

A(n - 1) ^ Vdf(i s ,I s Z) = 4(1 - n)(l + 5)d/ (Z) + 4(n + 1)C/(Z, V/), 

s=l 

which for n > 1 yields 

3 -, 

(3.9) 53 Vdf(H s ,I s Z) = -(1 + S)df(Z) + ^-—U(Z, V/). 

n — 1 

s— 1 

The sum of (3.7) and (3.8) gives 

3 

(3.10) (2n + 1) 53 V<tf(tsJaZ) = (2n + 1)(1 - 2S)df(Z) - (2n + 3)T°(Z, V/) - (2n + 5)C/(Z, V/). 

s=l 

Equalities (3.9) and (3.10) imply (3.5). Letting Z = J S V/ in the latter it follows T°(/ S V/,V/) = since 
C/(/ s V/, V/) = 0. □ 

3.2. Formulas for the derivatives of /. By assumption, the second order horizontal derivatives of / 
satisfy the Hessian equation (1.6). We derive next formulas for the second order derivatives involving a 
horizontal and a vertical directions. 

Lemma 3.3. With the assumptions of Theorem 1.3 we have 



2?i + 3 



(3.11) Vdf(^I 1 Z) = -df(Z) + —^—T (Z,Vf)-T (I l Z,hVf) + \ r MZ,Vf) 

4(2n+l)L J (2n+l)(n — 1) 



2n z + 3n - 1 



mi 



d 



n+1 



2n + l 



(3.12) Vdf(Z, eO = df(IiZ) - — I— T\hZ, V/) + T°(Z, /.V/) - - ■ —. -U(I t Z, V/) 



An 



(2n+l)(n - 1) 
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Proof. The second equality of (3.3) can be written in the form 



(3.13) AnQiihZ^f) = -df{Z) + {An-2)T{^,I i Z,^f) + Y d T{^,I s Z,Wf) 

s=l 
3 

+ (An + 2)Vdf&,IiZ) - ]T Vdf({ s ,I s Z) 



= -df(Z) + (An - 2) [-(T°(Z, V/) - T°(I t Z, /<V/)) - U(Z, V/)J + T°(Z, V/) - 3U(Z, V/) 



(1 + S)df(Z) 



n + 1 



U(Z,Vf) + (An + 2)Vdf(^hZ), 



n - 1 

where we used (2.9) and (3.9). Now, equalities (3.13), (3.5) and the third equality in (2.10) imply 

1 



(3.14) ^6,1,7) = -^-^ 



-df(Z) 



2n + 3 



t°(z,v/) + t°(/^,W) 

T°(Z, Vf)-T°(I t Z,I t \7f) 



(2n+l)(n - 1) 
2n 2 + 3n - 1 



4(2n+l)L v '' JJ K " % '->•■"] • (2n + l)(n-l) 
Finally, the Ricci identity, (2.9) and (3.11) yield 

(3.15) V 2 /(Z,&) =Vdf(£ i ,Z)+T(t i ,Z,Vf) 

= U(I t Z) + —p-^——T Q> (I i Z, V/) - p±T\Z,hVf) + J n l~u~\ mz,Vf) 
2 2(2n + 1) 2n + 1 (2n+l)(n— 1) 



U(Z,\7f) 
U(Z,Vf). 



df(hZ) 



n + 1 
2n+ 1 



T°(I t Z, v/) + r (z,/,v/) 



4/? 



(2n+ l)(n - 1) 



U(hZ,Vf), 



which completes the proof. □ 

Next, we compute the second vertical derivatives of /. We start with a basic useful identity involving 
only vertical derivatives. 

Lemma 3.4. With the assumptions of Theorem 1.3 the following identity holds 

n + 1 



(3.16) V/(6,6) = -/- 



An(2n + 1) 



(V ea T°)(e a ,V/) - (V ea T )(he a ,hVf) 



1 



(2n + l)(n- 1) 



(V eo C0(e„,V/). 



Proof. Differentiating (3.12), using (1.6) and (2.3) we obtain 

(3.17) V 3 f(X 7 Y, 6) - OLi (X) V 2 f(Y, 60 + a,(X)V 2 /(F, fc) 
n+1 
~2n + 1 



(VxT°)(/ i y, v/) + (v^T°)(y,/<v/) 

n+1 



4?? 



■(VxU)(IiZ,Vf) 



(2n + l)(n-l) 

(x,^y) + r°(/^,y)j+ 4n — -t/pr^y)} 

J (2vi + lj(n — 1) J 

+ dm,) { - ff (x, f) + — [to ( /^, / 4 y) - nx, y)] + (2w+1)(n _ 1) ^, n } 



2n+ 1 

n + 1 r 



+ dm {^Y) + ^[n Wi Y)+ni k X,Y)} + *" U(X,I>Y)} 



2n+ 1 
n + 1 



+ dm) { - Uj (X,Y) + ^ [T°( Wi Y) - T°( W) ] - (2 „ + i ;; n _ 1) ^^)} 

-U(hYVf) 



aj (X)\df(I k Y) - £±±T°(I k Y,Vf) - £±\T°(Y,I k Vf) 



2n+l 



2n+ 1 



OLk 



n+1 



n + 1 



(X)[df(I,Y) ^—T^L^f) ——^FJFf) 



(2n + l)(n - 1) 
4n 

(2n+ 7)(n - 1)' 
4n 

(2n +lj(n-l)' 
An 

(2n+ l)(n - 1)' 
4n 



(2n+ l)(n - 1) 



U(I 3 Y,Vf) 
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Applying again (3.12) to the second and the third terms in the first line we see that the terms involving the 
connection 1-forms cancel and (3.17) takes the following form 



(3.18) V 3 f(X,Y,^) 



n+1 



2n + 1 



(V A -T°)(/ 2 y, V/) + (\7 x T°)(Y,I t \7f) 
n + 1 



An 



2n + 1 

n + 1 r 



+m){- 9 {x,Y). 2n+1 



T°(x,/ l y) + T°(/ i x,y) 

T a (I t X, I{Y) - T°(X, Y) 



(2n + l)(n - 1) 
An 



(V x U){hZ,Vf) 
-U(XJiY)} 



2n+ 1 



T^IjXJ^+T^hX, Y) 
n + 1 



2n + l 



T^hXJ^-T^IjX, Y) 



(2n+l)(n - 1) 

(2n+l)(n-l) U{X ' Y) } 

— — — -U(X,I k Y)\ 

(2n + l)(n — 1) > 

An 



(2n+ l)(n - 1) 



-UiXJjY)}. 



On the other hand, the skew-symmetric part of (3.18) and the Ricci identity listed in the fourth line of (2.13) 
yield 

(3.19) V 3 f(X,Y,^)- V 3 f(Y,X,&) 
n+1 



2n+ 1 
An 



(VxT )(J,y, V/) + (VxT°)(r, JiV/) - (V Y T°)(I t X, V/) - (V Y T°)(X,I t Vf) 

An 



(2n + l)(n-l) 

+ 2df(^){u Jk (X,Y) 



(y x u)(i t Y,vf) - (v Y u){iiX,vf) 

n + 1 r 



2/ 



w<(A-,y) + 



+ 2d/(a){-o; j (X ! y) + 



2n+ 1 
n + 1 [ 



2n+l 



T°(I k X,Y)-T°(X,I k Y) 
T°(X,IjY) -T°{I 3 X,Y) 



(2n+ l)(n - 1) 



(2n+l)(n- 1) 
4n 



(2n + l)(n - 1) 



-c/(x,/,y)} 



-2df(^)p k (X, Y) + 2df(t k ) Pj (X, Y)-2^2 oj s (X, y)V 2 /fe, &). 



The trace X = e a ,Y = lie a of (3.19) and the second equality of (2.10) give (3.16), which completes the 
proof. □ 

Remark 3.5. The detailed proof of (3.18) shows a particular consequence of (2.3) which is that a covariant 
derivative of identities that are not Sp(l) invariant can lead to formulas which do not involve the connection 
one-forms. In the rest of the paper we shall usually skip many straightforward calculations some of which 
rely on a similar use of (2.3). 

3.3. The elliptic eigenvalue problem. In this sub-section we will show that (1.6) implies that / satisfies 
an elliptic PDE. Let A h be the Riemannian Laplacian of the metric (1.5). 

Lemma 3.6. On a qc manifold of dimension bigger than seven any smooth function satisfying (1.6) obeys 
the following identity 

n+1 ,„ „ nw . 3 



(3.20) 



A h f = (An + 3)/ + 



(V eo T )(e o ,V/) + 



: (V ea U)(e a ,Vf). 



n(2n+l) (2n+l)(n-l) 

Proof. It is shown in [40, Lemma 5.1] that the Riemannian Laplacian A h and the sub-Laplacian A of a 
smooth function / are connected by 

3 

(3.21) A h f = A/-^V 2 /(£ S ,£ S ). 

s=l 

Equation (3.21) is a consequence of the formula (2.2), A h f = - E«=i V h d/(e ,e )-£* =1 V h #(6,&), and 
the identities T(e a ,A 7 e a ) — T(£ s , A,£ s ) = 0, A £ T(TM) which follow from the properties of the torsion 
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(3.22) 



n + 1 
n{2n+ 1) 



(V e T°)(e a ,Vf) 



(2n+ l)(n - 1) 



(V ea C/)(e a ,V/). 



A substitution of (3.22) in (3.21), taking into account that / satisfies (1.6) hence A/ = 4nf, we obtain 
(3.20) which proves the lemma. □ 

A consequence of Lemma 3.6 and Aronsajn's unique continuation result, [2], is that |V/| cannot vanish 
on any open set. We note this important fact in the next remark. 

Remark 3.7. If M and f are as in Theorem 1.3 then |V/| ^ in a dense set since f ^ const. 

3.4. Formulas for the torsion tensors. In this sub-section we derive formulas for the components T° and 
U of the torsion tensor. 

Lemma 3.8. With the assumption of Theorem 1.3 the following identities hold true for any X,Y, Z 6 H 
(3.23) T°(7 s V/,/ t V/) = 0, s^t, s,t €{1,2,3}, 



(3.24) T°(V/,V/) 



6n 

n-1 



£/(V/,V/), T°(J S V/,J S V/) 



2n 
n - 1 



f/(V/,V/), se {1,2,3}, 



(3.25) 



V/| 2 T U (Z,V/) 



6n 
n- 1 



J7(V/, Vf)df(Z), |V/| 2 [/(Z, V/) = £/(V/, Vf)df(Z), 



(3.26) 



2r? 

|V/| 4 T (X,F) = -C/(V/, V/) idf{X)df{Y) -Y^df{I s X)df{I s Y) 



1 3 

(3.27) |V/| 4 C/(Z,X) = -{/(V/, V/) [iV/pfffZ,!) - + V rf/(I s Z)rf/(/ s l) 

n — 1 L V *— ' 

s— 1 

Proof. To determine the torsion tensors T° and {/ we are going to apply the following identity [42, 44] for 
the [—1] component of the curvature 

(3.28) 3R(Z, X, Y, V/) - R(hZ, I X X, Y, V/) - R(I 2 Z, I 2 X, Y, V/) - R(I 3 Z, I 3 X, Y, V/) 

= 2 [g(X, Y)T°(Z, V/) + g(Z, Wf)T°(Y, X) - g(Y, Z)T°(X, V/) - 5 (V/, X)T°(Y, Z) 

3 

-2^ [u s (X,Y)T°(Z,I s Vf)+uj s (Z, Vf)T\X,I s Y) -u t (Z,Y)T°(X,I,Vf) -u s (X,Vf)T°(Z,I s Y) 

s=l 

3 

+ [2u} s (Z,X)(T°(Y,I s Wf] - T°(I S Y, V/)) - 8w s (F, Wf)U(I s Z,X) - 4Su s (Z,X)u s (Y,Vf) . 

s=l 

With the help of the Ricci identity, cf. the second equality of (2.13), we write the curvature tensor given by 
(3.2) in the form 



(3.29) R(Z,X,Y,Vf)= df(Z)g(X,Y)-df(X)g(Z,Y) 



[vdf(Z, £ s )u s (X, Y) - Vdf(X, Z s )lo s (Z, Y) - 2Vdf(Z s ,Y)w s (Z, X) 



s=l 
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A calculation shows 
3 



t=i 



(3.30) J2R(I t Z,I t X,Y, V/) = \df{I s Z)uj s {X,Y) - df(I s X)uj s (Z, Y) 

s=l 

3 

J2 \vdf{I t Z^ s )uj s {I t XX) - Vdf(I t X,£ s )cj s (I t Z,Y) - 2Vdf(£ s ,Y)uj s (I t ZJ t X) 

s,t=l 

3 

= [df(I s Z)u> s {X, Y) - df(I s X)uj s (Z, Y) + 2Vdf(t s ,Y)oj s (Z, X) 



s=l 



- g(X, Y) Vdf(I s Z, £,) + g(Z, Y) £ Vdf(I s X, £„) - £ Ui (X, Y) [vdf^Z, &) - Vd/(/ fc Z, 

S=l 8=1 (ijfe) 

+ ^ Wi(Z, Y) [vdf(I 3 X, &) - Vdf(I k X, &)] , 

where X)(yfc) denotes the cyclic sum. Now. (3.29) and (3.30) together with (3.11) and (3.12) yield 
(3.31) 3R(Z, X, Y, V/) - R(hZ, I,X, Y, V/) - R(I 2 Z, I 2 X, Y, V/) - R(I 3 Z, I 3 X, Y, V/) 

3 3 3 

= g(X, Y) hdf(Z) + Y V 2 f(I s Z, &)] - 9(Z, Y) hdf(X) + Y V 2 f(I s X, - 8 J] lo s (Z, X)Vdf{^ sl Y) 



s=l 



8 = 1 



]T Ui(X,Y) hv 2 f(Z^i) - df{hZ) + V 2 /(-^,60 - V 2 f(I k Z,^) 



(ijk) 



Y "i(Z,Y) [3V 2 f(X,^) - df(UX) + V 2 f(IjX,Z k ) - \7 2 f(I k X,^ 

(ijk) 



= g(x,Y) 
-g(z,Y) 

3 

Y^s{Z,X)Usdf{I s Y) + 



An + 4 



T°(Z,Vf) 



Yin 



-U(Z,Vf) 



2n + l~ v ' " J ' ' (2n + l)(n - 1) 
in + 6 



2n+ 1 

' 4n + 4 



(2n+l)(n - 1) 
T°(I s Y,Vf)-T Q (Y,I s Vf) 



(2n+ l)(n- 1) 

E N r 4n + 4 n , N 4n . , ~i 

Y i^Ti T ^ ^ + (2n + iKn - 1) ^ (/ ^' V/) ) 



U{I s Y,Vf) 



(2n + l)(n - 1) 
-4 « + 4 ^o. 



X>(^){^ r °(*>w) 



4n 



s=l 



-C/(7 S X,V/)}. 



(2n + l)(n - 1) 

Subtracting (3.28) from (3.31) and applying (3.11), (3.12) and the properties of the torsion we come to 



(3.32) = g(X, Y) \t°(Z, V/) + -^-U(Z, V f )j - F) [r°(X, V/) 
L n — 1 J L 



n- 1 



U(X,Vf) 



-y>,(*,y) r°(z,/ s v/) + — -c/(/ s z,v/) + 5>.(z,y) t°(x,i.v/) + — -u(i s x,vf) 

* — ' L 71 — 1 J * — ' L 71—1 

s=l s=l 



- V oj 8 (Z, X) \2T°(I 8 Y, V/) - 2T°(Y, 7 8 V/) -U(I S Y, Vf) 

*— ' L n — 1 

s— 1 
3 

(2n + 1) Y [df(I s X)T°(Z, I S Y) - df(I s Z)T°(X, I S Y) - 4df(I s Y)U(I s Z, X) 



+ (2n + l)d/(X)T°(Z, y) - (2n + l)4f (Z)T°(*, F). 
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(3.33) (2n + l)\Vf\ 2 T°(X,Y) = (2n + l)df(X)T°(Wf, Y) 
+ 5 (X,F)[t°(V/,V/) 

3n 2 - 2n - 4 



' t/(V/,V/)l -df{Y)\T°(X,Vf) + J^-U(X,Vf) 



n — 1 

3 

^d/(/ s y)[T°(x,/ s v/) 



n-l 



n — 1 
U(I s X,Wf) 



- J2 d f(^ x ) \2T\Y, I s Vf) + (2n - 1)T°(/ S y, V/) 



s+l 



n- 1 



U{I a Y,Vf) 



Letting Y = V/ in (3.33), then using (3.34) and (3.6) shows 

In 

(3.34) |V/| 2 T°(X, V/) = T°(V/, V/)4f (*) ' 



(n+l)(n-l) 

On the other hand, letting X = I\V f in (3.33), using (3.6) and (3.34) gives 

8rc 2 - 8ra - 4 



U(Vf,Vf)df(X) - |V/| 2 C/(X, V/) 



(3.35) = -df(hY) T°(V/, V/) + T°(JiV/, /i V/) - 



n- 1 



f/(v/,v/) 



- d.f(I 2 Y)T (hVf,I 2 V.f) - df(hY)T°(hVfJ^f) 

- (2n - 1)|V/| 2 [t°(Y, /i V/) - T°(7iY V/)] + |v/f 

From (3.35) with Y = 7 2 V/ and (3.8) the identity (3.23) follows since |V/| 2 ^ 0. 
Setting Y = I\\7 f into (3.35) implies 



n- 1 



U(hY,S7f). 



(3.36) 



4/; 



T°(V/,V/) + T°(IiV/,JiV/) = T f/(V/,V/) 

n — 1 



The latter equality together with the (2.8) yield (3.24). 
The equalities (3.35), (3.23) and (3.24) imply 

(3.37) (2n - 1)|V/| 2 T°(Y,7 S V/) = (2n - 1)|V/| 2 T (7 S Y V/) + |V/ 



n - 1 



t/(7 s Y,V/) 



df(I s Y) [T°(V/, V/) + T°(7 S V/, J.V/) - 8 " V (V/, V/) 



n — 1 



(2n-l)|V/| 2 T°(J s F,V/) + |V/p 



n — 1 



C/(7 S Y, V/) + 4(2n + l)df(I s Y)U(Vf, V/). 



Let Y" = TiV/ in (3.32) in order to see 



(3.38) 2(2n+l)\Vf\ 2 U(hZ,X) = wi(Z,A-)[t°(V/, V/) +T°(7 1 V/,7 1 V/) - 

+ nd/(X)[r°(Z,7 1 V/) - -J—^Z.V/)! -nd/^fT^X^iV/) - 
L n — 1 J L 

+ n#(JiX) [t°(Z, V/) - -^—U(Z, V/)l - nd/(7xZ) [t°(X, V/) - 



ndf{I 2 X) T°(Z,I 3 Vf) 



n — 1 
1 



n — 1 



U(I 3 Z,Vf) -ndf{hZ) T\X,hV.f) 



n 




— 




1 




n 




— 




3 




n 




— 






1 









U(hX,Vf) 
U{X,Vf) 
-U{hX,Vf) 



- ndf(hX) \T (Z, 7 2 V/) - -!—U(I 2 Z, V/)j + nd/(7 3 Z) [t°(X, 7 2 V/) - -*—U{I 2 X, V/) 
L n — 1 J L n — 1 

Letting X = V/ in (3.38) and applying (3.24) we obtain 

(3.39) (4n 2 - n - 2)U(Z, V/)|V/| 2 = (6n 2 - n - 2)C/(V/, Vf)df(Z) - n(n - l)T°(hZ, hV '/)|V \f\ 2 . 
Therefore, 



(3.40) 



|V/| 2 n(n - 1)T°(Z, V/) - 3(4n 2 - n - 2)U(Z, V/) = -3(6n 2 - n - 2)C/(V/, Vf)df(Z). 
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On the other hand, taking into account (3.24), equality (3.34) yields 

(3.41) | V/| 2 [(n 2 - 1)T°(Z, V/) + 3nU(Z, V/)] = -3n(2n + 1)£/(V/, Vf)df(Z). 

Solving the system of equations (3.40) and (3.41), we obtain (3.25). 

A substitution of (3.25) and (3.43) in (3.33) gives (3.26). Now, a substitution of (3.25) and (3.43) in 
(3.38) shows (3.27). □ 



We finish this section with a few useful facts. As a direct corollary from (3.5), (3.25) and (3.43) it follows 
(3.42) 

In addition, from (3.25) and (3.37) it follows 



| V /| 2 (5-2) = ^^f/(V/,V/). 



2n 2n 
(3.43) |V/| 2 T°(Z,/ S V/) = — -U(Vf,Vf)df(I s Z), |V f\ 2 T°(I s Y, / a V /) = -U(Vf,Vf)df(Z). 



n — 1 

The equalities (3.25) and (3.43) yield 



n — 1 



(3.44) T Q (I s Z,Wf) = 3T°(Z,I s Vf), T°(Z,Vf) = -3T°(I s ZJ s Vf), T°(Z,V/) 



6n 
n — 1 



3.5. Formulas for the covariant derivatives of the torsion tensors. Here we shall prove formulas for 
the covariant derivative of the torsion tensor. 

Lemma 3.9. If M and f are as in Theorem 1.3, then we have the following identities for the covariant 
derivatives of the torsion tensor at the points where |V/| ^ 0, 

(3.45) |V/| 2 (V Z T°)(A, Y) = ^±1 fdf(Z)T°(X, Y) 



n + 2 



2n t/(V/,V/) 



n-1 |V/| : 



/ [ - 3df(Y)g(X, Z) - 3df(X)g(Y, Z) + £ (df(I a Y)u s (X, Z) + df(I s X)u s {Y, Z)) 



2n U(Vf,Vf) 
n-1 |V/| 2 



Udf(Y)u>i{X, Z) + df(IiY)g(X, Z) - df{I 3 Y)u k {X, Z) + df{I k Y)u 3 (X, Z) 



(ijk) 

2n t/(V/,V/) ^ ^ df{x)uJi{Yj z) + df ( IiX )g(Y, Z) - df{I 3 X)u k {Y, Z) + df{I k X)uj 3 (Y, Z) 

(ijk) 



-i |v/| 5 



(3.46) \Vf\ 2 (V z U)(X,Y) -2fdf(Z)U(X,Y) -2j2df^ s )df(I s Z)U(X,Y) = 



2n - 2 
n + 2 



■fdf(Z)U(X,Y) 



2j2df(ts)df(I s Z)U(X,Y)- 



1 C/(V/,V/) 



s=l 



n-1 |V/| 5 



- 2fdf(Z) -2j2df(t s )df(I s Z)]g(X,Y) 



8 = 1 



C/(V/,V/) 



n-1 |V/| 2 
n U(Vf,Vf) 



f[df(Y)g(X, Z) + df(X)g(Y, Z) + £ (df(I s Y)uj s (X, Z) + df(I s X)u s {Y, Z)) 



n-1 |V/| 2 

n t/(V/,V/) 
n - 1 |V/| 2 



J2 dfiti) k(y)^i(X, Z) - df(IiY)g(X, Z) + dfiljY^kiX, Z) - df(I k Y)cj 3 (X, Z) 



(ijk) 



^2 df{ii) \df{X)^{Y, Z) - df(IiX)g(Y, Z) + df(I 3 X)u; k {Y, Z) - df(I k X)cj 3 (Y, Z) 



(ijk) 



where X^yfc) means the cyclic sum, cf. Convention 1.5. 
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Proof. The contracted Bianchi identity reads [37, 44] 

(3.47) (V ea T°)(e a ,X) + ^±i(V etl C/)(e , X) - (2n + l)dS(X) = 0. 

n — 1 

After taking the trace in the covariant derivatives of (3.25) and (3.42) we obtain 



(3.48) 



(V e „E0(e fl 



o,v/) = V/ — ——5 — -4n/ — — -rTg — , 



4n + 4 



|V/| 
V/ 



c/(v/,v/) 



|v/| 



|v/| 



V/(5) 

The system (3.48) and (3.47) imply 

(349) n /(V/.V/) )-i"-l ^(V/.V/) 

ld 9j V/ V |V/| 2 J n + 2 ; |V/| 2 

Similarly, using in addition (3.44), we have 
(3.50) 



(V eo C/)(e ,7 s V/) =I S V/^ — p^yp — ) + 4nd/(£ s ' 



|V/| 





n — 1 



|V/| 5 



|v/p 



/.V/(5) 



4n + 4 
n-1 



|V/|- 



Since the differentiation of (3.43) involves covariant derivatives of the almost complex structures the deriva- 
tion of (3.50) requires some care we do it explicitly again, cf. Remark 3.5. We start with the proof of the 
first formula in (3.50). Differentiating the first equation in (3.43), taking into account (2.3), we have 

(V X T°) (Z,/iV/) - aj (X)T°(Z,I k Vf) + a k {X)T°(Z,I j Vf) + T°(Z,I l V x (V/)) 



X 



2n 
n- 1 



2n C/(V/,V/) 



f U(Vf,S7f) 
\ |V/| 2 



df{hZ) + U[ ^J ] 9 (Vx (V/) , hZ) 



aj {X)df (I k Z) + a k (X)df (IjZ)] 



n-1 |V/| 2 
The formula for the Hessian (1.6) gives 

3 

(V X T°) (Z, IiVf) - aj {X)T°(Z, h V/) + a k (X)T°(Z, 7,-V/) - /T (/ 2 X, Z) - ^ # (&) T°(V«*, 



2n 



( U(Vf,Vf) 

V |v/| 2 



27i U(Vf,Vf) 
n-1 |V/| 2 



- aj (X)d/ (/ fc Z) + a*(X)df (7,-Z)] . 



Taking the trace in the above identity and then applying the first equation in (3.43) to the obtained equality 
we see that the terms involving the connection 1-forms cancel, which gives the first identity in (3.50). 

The second line in (3.50) follows similarly. 

The system (3.50) and (3.47) yields 



(3.51) 
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We calculate the divergence of T° differentiating (3.26), taking the trace in the obtained equality and applying 
(3.49), (3.51). After a short computation we obtain 



(3.52) |V/| 2 (V ea T°)(e a ,y)-2/T (r,V/) + 2^d/(e s )T (r,/ s V/) = 

s=l 

f { |v/| 2 ) /( } + 4j * v/ V |V/| 2 J /( s } 

3 

3V 2 /(e a , e a )df(Y) - £ V 2 /(e a , I s e a )df(I s Y) 

s=l 

2n C7(V/,V/) 



n — 1 

2n [/(V/,V/) 
n-1 |V/| 2 



n-1 |V/| 2 



3V 2 /(v/, y) + £ v 2 /(/ s v/, / s y) 



s=l 



Applying (1.6), (3.25), (3.43), (3.49) and (3.51) to (3.52), we get 



(3.53) |V/| 2 (V et T°)(e a ,y) 
12n ,tf(V/,V/) 



-1 J |V/| 2 



-i |v/| 2 



«=1 



n + 2 |V/| 2 



6n E/(V/,V/) ,^ 6 ? i Uiyf^Tf)* 
+ n-1 |V/| 2 - n-1 |V/| 2 V ' 

r 3 

- J^^T^ £ [-M(y) + df^)df(hY) - dfi&df&Y) - df(Z k )df(I k Y) 

12n(n + l)(2n+l) l/(V/,V/) 2n + 1 U(Vf,Vf) ^ 

= (n + 2)(n-l) iVfiS- M(r) " 4 " n - 1 iVfl 2 L^'WW 



Applying (3.25) and (3.43) to (3.53), we derive 



(3.54) |V/| 2 (V ea T°)( ea ,y) = - (4 " + ^ + 1} /T°(y, V/) + (4n + 2) V df(L)T a (Y, I.Vf). 

n + 2 ^— f 
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Now, we calculate the divergence of U differentiating (3.27), taking the trace in the obtained equality and 
applying (1.6), (3.49) and (3.51). We have 



(3.55) |V/| 2 (V ea lO(e Q , Y) - 2fU(Y, V/) + 2 £ df{£ s )U(Y, J S V/) 



1 „fU(Vf,Vf) 



-Y 



- r v |v/| 2 

1 C/(V/,V/) 



iv/p 



n - 1 



V/ 



|V/| 2 



df(Y)-Y,l^f 



n-1 |V/| : 



s=l 
3 



(V/, V/) 
|V/| 2 



2V 2 /(Y, V/) - nV 2 /(e a , e a )d/ (Y) - V 2 /(e a , I s e a )df(I s Y) 



n C/(V/,V/) 



-i |V/| 2 



V 2 /(V/,Y)-^V 2 /(/ s V/,/ s Y) 



1 v /f/(v/,v/) 



n — 1 



-Y 



2// 



|V/| 



2/? 



|V/| 



2n 
n — 1 



n - 1 



n + 2 



fU(Y, V/) - df(£ s )U(I s Y, V/) 



4n 2 - 2 



-/C/(Y, V/) + df($ s )U(I s Y, V/) /[/(Y, V/) + ^ df(Z s )U(I s Y, V/) 

1 rt — 1 : f n — 1 L ' 



Thus, from (3.55) we obtain 



(3.56) |V / | 2 (V ea t/)(e a ,Y) = --l T Y( ^ / /| y /) )|V/' 2 



4 " 2 + 6 " /C /(Y, V/) - ^-^ ± dm)U(I s Y, V/). 

71—1 



n + 2 



A substitution of (3.54), (3.56) and 



2wm _ 4n + 4 / £/(V/,V/ )^ |V7J ,, 2 
|V/| 2 



|V/| 2 Y(S) 



71 - 1 



|V/| 5 



in (3.47) implies 



(3.57) 



2n-2 



|V/| 2 



s=l 



The equalities (3.56) and (3.57) yield 



(3.58) 



|V/| 2 (V ea [/)(e a ,Y) = - 



2(n + l)(2n+ 1) 
n + 2 



/C/(Y V/) - 2(2n + 1) ^ df(QU(I s Y, V/). 



20 



S. IVANOV, A. PETKOV, AND D. VASSILEV 



We calculate from (3.26) using (1.6), (3.25) and (3.57) that 

3 

\Vf\ 2 (V z T°)(X, Y) - 2fdf(Z)T°(X > Y)-2J2 df(£ s )df(I s Z)T°(X, Y) 

2n - 2 



s=l 



fdf(Z)T°(X,Y) -2Y^df{Qdf{I s Z)T\X,Y) 

s—1 

2n t/(V/,V/) 



n-l |V/| 2 

2n t/(Vf,V/) 
n-l |Vf| 2 



2n t/(Vf,V/) 
n - 1 |V/[ 2 

2n l/(V/,V/) 
n-l |V/| 2 

2n £/(V/,V/) 
n-l |V/| 2 

2n C/(V/,V/) 
n-l |V/| 2 

2n UjyjjVf) 
n-l |V/| 2 
2n C/(V/,V/) 



/ [ - 3d/(y) 5 (X, Z) + £ Z) 

s=l 
3 

/ [ - 3df (X) 5 (F, Z) + ^ # (Z.X)w.(Y; Z) 

s=l 

#(6) [3# (*>i(*, Z) + 4f(7 1 y) fl (X, Z) - df (7 2 Y)a;3(X J Z) + d/(J 3 y)a; 2 (X, Z) 
[3tjf(A>i(y, Z) + df(hX)g(Y, Z) - d/(/ 2 J5£> 3 (y Z) + df (I 3 X)uj 2 (Y, Z) 
df (6) [3df (y)w 2 (X, Z) + df(I 2 Y)g(X, Z) + df(hY)u 3 (X, Z) - d/(J 3 y)o;i(X, Z) 
4f (fa) [3d/(X)o; 2 (y, Z) + df(I 2 X)g(Y, Z) + df(hX)uj 3 (Y, Z) - df (hX^Y, Z) 



4f (&) [34f(y)w3(x, z) + df(i 3 Y) g (x, z) - d/t^y^cx, z) + # (j 2 i>i(*, z) 

[34f (X)w 3 (y, Z) + df(I 3 X)g(Y, Z) - rf/(/ 1 X) W2 (y Z) + df{I 2 X)w x {Y, Z) 



n-l |V/| 2 

The last equality yields (3.45). Finally, equation (3.46) follows from (3.27) using (1.6), (3.25) and (3.57). 

In the next, key step of the proof, where we show that the torsion tensor vanishes, we shall use 
following particular cases of Lemma 3.9. For Z = V/, (3.45) gives 

2n — 2 r i 

(3.59) (v V /T°)(x,y) = __ / T (x,y) + df(c 1 )[T (/ 1 x,y) + r (x,/ 1 y) 

+ df{&) [T (I 2 X, Y) + T°(X, 7 2 y)] + # (&) [T (I 3 X, Y) + T°(X, I 3 Y) 
Similarly, letting Z = 7 V/ in (3.45) we obtain 

(3.60) (V 7iV fT°)(X, Y) = 2df(£ i )T°(X, Y) + f \T Q {I t X, Y) + T°(X, hY) 



df(Q T°(I k X,Y) + T°(X,I k Y) 



T a {IjX, Y) + T Q (X,IjY) 



The substitution of Y = V/ in (3.54) taking into account(3.24) and Lemma 3.2 gives 
(3.61) 

while the substitution Z = e a , X = 7e a , Y = ijV/ in (3.45) and (3.43) gives 
(3.62) 



|V/r ( V eo T° )(eaj V/) = ^(n + l)(2n + l) 

(n + 2)(n - 1) 



V/| a (Ve.T°)(W 4 V/) = _ 4n(n + l)(2n + l) /[/(v/;V/) ^ 



(n + 2) (n-l) 

Finally, letting Z = V/, I s Vf in (3.46) shows the next equality 

2n - 2 



(3.63) 



(y vf U)(X,Y) 



n + 2 



fU(X, Y), (V Is v f U)(X, Y) = 2df(^ s )U(X, Y). 



THE OBATA SPHERE THEOREMS ON A QUATERNIONIC CONTACT MANIFOLD 



21 



3.6. Vanishing of the torsion. In this section we show the vanishing of the torsion, T° = U = 0, by 
calculating in two ways the mixed third covariant derivatives of a function satisfying (1.6). 

Lemma 3.10. If M satisfies the assumptions of Theorem 1.3, then the torsion tensor vanishes, T° = 0, 
U = 0, i.e., M is a qc-Einstein manifold. 

The proof occupies the remaining part of this sub-section. 



3.6.1. The Ricci identities. We are going to use the sixth line in (2.13). A substitution of the contracted 
Bianchi identity (3.47) in the second formula of (2.12) gives 

(3.64) (2n + l)p i (^,I k X) = -(2n+ l)ft(&>^*) 

1 n 

= — A [^eT ){e a ,X) + {V ea T ){Iie a ,IiX)} + ~{V ea U)(e a ,X). 

4 n—1 

Let Z = V/ in (2.11), and then substitute the obtained equality in the sixth formula of (2.13), after which 
use (3.64) in order to see 



(3.65) V 3 /te,X,Y)-V 3 /(X,y,&) 

= -V 2 /(T(6, X), Y) - V 2 f(X, T(&, Y)) - df((\7 x T)(^ Y j) + (V x U)(IiY, V/) 

+ \ [(VyT°)(7V/,X) + (VyT°)(V/, JiJC)] - \ [(V V /T°)(7y,X) + (V V /T°)(y, 



1 



2n+ 1 
1 

2n + 1 



1 n 

■-(V eo T°) [(e ,I fe V/) - (I k e a ,Vf)} + -(V e „[/)(e fl ,4V/) 

4 n—1 



1 



1 



2n4 


1 


1 




2n4 


1 


1 




2n + 


1 


1 




2n + 


1 



{W ea T°)((e a ,I k Y)-(I k e a ,Y)) 



-(V e T°)«e ai I 3 Y) - (I,e a ,Y)) - -—(V ea U){e a ,I Y) 



(X7 ea T°){{e a ,I k X)-(I k e a ,X)) 



n - 


- 1 


n 




n — 


T' 


n 




n — 


T< 


n 




n — 


i« 


n 


t 



{V ea U){e a ,I k Y) 



-(V ea T°) «e ai I 3 X) - (I 3 e a ,X)) - — -(V ea U)(e a , IjX) 



Uj{X,Y) 
u k {X,Y) 
df(IjX) 
df(hX) 
df(IjY) 
df(hY). 



Note that from (2.9) we have 
(3.66) 



T(Zi,X,Y) = -1 [T a (I t X,Y) + T°(X,I t Y)] -U(X,hY). 



Differentiating the above formula we find, applying (2.3), 



(3.67) df((V x T)(^,Y)) = -- [(VxT )(I i y ) V/) + (V JC T )(y,J i V/)] + {V X U){UY,V f). 
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Using (1.6), (3.66), (3.67) and the properties of torsion tensor listed in (2.8), we obtain from (3.65) 
(3.68) V 3 f(ti,X,Y)-V 3 f(X,Y,Zi) 

= \ [(VxT°)(IiY,Vf) + (VxT°)(y,X t V/)] + \ [(V Y T°)(IiX,Vf) + (V Y T°)(X, I<V/)] 
- \ [(Vv/T°)(X, hY) + (V V /T°)(7iX, y)] - i [T^JiX, y) + T°(X, IiY)] f 



i (t°(x, i fe y) - T°(/ fe x, y)) + 2t/(x, j fc y) 

i (r°(/,X, y) - T°(X, I 3 Y)) + 2U(I j X, Y) 



rf/fe) 



1 



2n+ 1 

1 

2n + 1 



- i (Ve £ ,T )[(e a ,/ J V/)-(/ J e a ,V/)] 



1 



2n + 


1 


1 




2n4 


1 


1 




2n + 


1 


1 





^(Ve„T°)((e a ,I fc y) - (he a ,Y)) - -^—(V ea U)(e a ,I k Y) 
{V e T»){(e a J 3 Y) - (I jea ,Y)) - -—{V ea U){e a J 3 Y) 



2n+l 



i(Ve o T ) ((e a ,I k X) - (I k e a ,X)) 
1 



(V ea T°) ((e^/.-X) - (i^X)) - -JL_(W ea U){e a ,IjX) 



n - 


1 


n 




n - 


1 


n 






I 1 


n — 


n 






T< 


n — 


n 




n — 


T< 


n 


t 



(W ea U)(e a ,I k Wf) 
(V ea [/)(e a ,/ J V/) 



(V ea U){e a ,I k X) 



Uj{X,Y) 
to k {X,Y) 
df(IjX) 
df(hX) 
df(IjY) 
df(hY). 



For X = IiVf, Y = V/ equation (3.68) together with (3.59), (3.60), (3.6) and (3.23) imply 

(3.69) v-'y.-./.v/. V/) - v 3 /(/ 4 v/, V/,6) = o. 

On the other hand, a subtraction of (3.18) from (3.1) with A = £j gives 

71+1 



(3.70) V 3 f^ h X,Y)-W 3 f(X,Y,^) 
4n 



2n+ 1 
?? 4- 1 

— — — — -(v Y c/)(/,y, v/) - f-±— 

(2n + l)(n — 1) 2n + 1 

Uo 



(v x r°)(/y, v/) + (v x T°)(y, i&f) 

i An 

r (x,/ l y) + T°(/ i x,y) 



- dftti) 

dm) 



2n + l 

71 + 1 

2n + 1 

71 + 1 



271+1 



T°(IiX,IiY)-T (X, Y) 

r°(/ J x,/ J y) + T°(/ fe x, y) 
T°(/ fc x,/iy) -T°(/jX,y) 



4n 


(2nH 


-l)(n- 


1) 




4ti 




(2ra- 


-l)(n- 


1) 




4n 




(2n- 


hl)(n- 


1) 



(2n + l)(n- 1) 
dm)U(X,Y) 

df(^)U(X,I k Y) 

dfi^UiXJjY) 



fU{X,kY) 



£[V a /(&,&) + /] u.(X,Y). 



Letting X = ijV/, y = V/ in (3.70) and then applying (3.6) and (3.23) we obtain 

(3.71) V 3 /(&, W,V/)-V 3 /(W,V/,&) - ^"tlV /.v/^XW, V/) 

2n + 1 

4n 



(2n + l)(n - 1) 



(V/ iV /[/)(^V/,V/)-|^-/ 

4n 



(27i+l)(n - 1) 



T°(/ 4 V/,/ 4 V/)-T°(V/,V/) 

/t/(V/,V/)+[V 2 /(&,£i) + /] |V/| 2 . 
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Using (3.60), (3.63) as well as (3.24) in (3.71) we conclude 

(3.72) v 3 /(&, w, V/) - v 3 /( W, v/, eo = — -fuiyf, V/) + [v 2 /(6, 6) + /] I v/| 2 . 

n — 1 

The formula for the last term is given in (3.16) to whose right-hand side we apply (3.61), (3.62) and (3.58) 
in order to obtain 

(3 73) V 2 f (£■ £■) + f - 2(n+l)(2n+l) Cf(V/,V/) 

^■ 7dj V/lW,,+; " (n + 2)(n-l) ; |V/| 2 

Now (3.73) applied to (3.72) allows us to conclude 

(3.74) V 3 /(6,iiV/,V/)- V 3 /(JiV/,V/,fc) = _^/tf(V/,V/). 

n + 2 

Comparing (3.69) and (3.74) we obtain fU(\7f, V/) = 0, which implies U(\7f, V/) = taking into account 
Remark 3.7. Hence, T° = {/ = due to (3.26) and (3.27). This completes the proof of Lemma 3.10. 

3.7. The Riemannian Hessian. Here we show that if T° = U = the equality (1.6) implies that the 
Ricmannian Hessian satisfies (1.1) and therefore the manifold is the standard sphere due to the Obata's 
theorem. 

Lemma 3.11. Let (M, r), g, Q) be a qc-Einstein manifold, T° = U = 0, of dimension An + 3 > 7. Let h be 
the associated Riemannian metric (1.5). If f is smooth function whose horizontal Hessian satisfies (1.6), 
then the Riemannian Hessian of f with respect to the metric h satisfies (1.1). 

Proof. Taking into account (2.2) we have the following formula relating the Hessian with respect to the 
Levi-Civita and the Biquard connections 

(3.75) (V h ) 2 f(A, B) = V 2 f(A, B) + ± [h(T(A, B), df) - h(T(B, df),A) + h(T(df, A), B)] ,A,B e T(TM). 
From (3.75), (2.5) and (1.6) it follows that 

(3.76) {V h ) 2 f(X,Y) = -fh(X,Y). 

Let us recall that a qc-Einstein manifold, T° — U — 0, has integrable vertical space [37] thus the fourth line 
in (2.10) shows 

(3.77) h(T(Z.,Z t ),X) = 0, 
Now, using (3.12) with T° = U = 0, we calculate from (3.75) 



(3.78) (V' l ) 2 /(X,&) = df{hX) + i/i(T(X,&), V/) - ±h(T{ti,Vf),X) - \h(T&, ]T X) 

S — l 

1 1 3 

+ -h(T(yf,X),& + -HT(J2df(L)L,X),Zi) = df{hX) +oj i (Vf,X) = 0, 

8=1 

taking into account (3.77) and the properties of the torsion (2.9) and (2.5). A similar argument shows the 
identity 

(3.79) (V fe ) 2 /(&,&) = V 2 /(6,&) = -/, 

where we have used (3.16) taken with T° = U = 0. 

Finally, we have to show (V /l ) 2 /(£i, £j) = 0. The trace with respect to X = e a , Y = Ije a in (3.19) together 
with the second equality in (2.10) and the condition T° = U = yields 

(3.80) vVfe,ei) = (i-sM/(a). 

Now, the equality (3.75) together with the fourth equality in (2.10), (3.77) and (3.80) imply 

(3.81) (V l ) 2 /te,G) = (i - s)4f(&) + is4f(G0 = a - \s)df^ k ) = 0, 

since (3.5) shows S = 2 in the case T° = U = 0. □ 
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At this point, applying the Obata theorem we conclude that our manifold is isometric to the unit sphere. 
In order to show that it is qc-equivalent to the sphere we shall use a Liouville-type result in the quaternionic 
contact case which we present next. 

3.8. Quaternionic contact conformally flat manifolds and the Liouville theorem. We start by 
recalling some basic definitions and facts. 

Definition 3.12 ([37, Section 7.2]). A diffeomorphism F of a qc manifold (M,rj) is called a conformal 
quaternionic contact automorphism (abbr. conformal qc- automorphism or qc-conformal map) if F preserves 
the qc structure, i.e. F*rj = [ity • 77, for some positive smooth function (i and matrix $ G 50(3) with smooth 
functions as entries, where r/ = 7?2) %)' is a local 1-form considered as an element 0/R 3 . If fi = const 
we call F a quaternionic contact homothety (abbr. qc-homothety). Finally, if /! = 1 then F is called a 
quaternionic contact automorphism (abbr. qc- automorphism). 

The above definition extends in the obvious manner to a map between two qc manifolds giving the notion 
of a qc-conformal map. The qc-conformal curvature tensor W qc , introduced in [42], is the obstruction for 
a qc structure to be locally qc-conformal to the fiat structure on the quaternionic Hciscnbcrg group G (H) . 
A qc-conformally flat structure is also locally qc-conformal to the standard 3-Sasaki sphere due to the local 
qc-conformal equivalence of the standard 3-Sasakian structure on the (4n + 3)-dimensional sphere and the 
quaternionic Heisenberg group [37, 42] via the Cayley transform whose definition we recall below, see also 
[44, Section 2.3.1] for the history and references of the Cayley transform on groups of Heisenberg type. In 
view of the uniqueness of the possible associated almost complex structures [6], (see also [37, Lemma 2.2]), 
a conformal quaternionic contact automorphism will also preserve, unlike the situation in the CR case, the 
associated almost complex structures and the conformal class [g] of the horizontal metric g on H . 

We shall use the following model of the quaternionic Heisenberg group G (H), see [37, Section 5.2] or [44, 
Section 4.3.4]. Define G (H) = i" x Iml with the group law given by 

(3.82) (q ,u> ) o (q,u>) = (q Q + q,uj + uj + 2Imq q), 

where q, q Q G H™, oj,u G ImH and H is the space of quaternions. We shall denote with (q,uS) the elements 
of G (H). The standard 3-contact form — (©1, ©2, 63) written as a purely imaginary quaternion valued 
1-form is defined by 

(3.83) 9 = i {duo - q ■ dq + dq ■ q) . 

Let 5 4 "+ 3 = {|g'| 2 + \p'\ 2 = 1} C H n x H be the unit sphere in the (n + l)-dimensional quaternion space 
equipped with the qc structure defined by the standard contact form fj on the sphere, 

(3.84) fj = dq' ■ q' + dp' ■ p' - q' ■ dq 1 - p' ■ dp'. 

Note that fj is twice the standard 3-Sasakain form on the sphere. Following [37], see also [44], we identify 
the quaternionic Heiseneberg group G (H) with the boundary of a Siegel domain in IP x H, 

G (H) = {yei"xi : D\t p = | 9 | 2 }, 

using the map (q,w) n> (q, \q\ 2 — lo). Under this identification we have 

(dp — dp) + 2dq ■ q — 2q ■ dq . 

The Cayley transform, see [52] and [21], is the map 6 : S^ 3 \ {(-1,0)} -> G (H), (-1,0) G H" x H, from 
the sphere S 4n+3 = {\q'\ 2 + \p'\ 2 = 1} C H" x H minus a point to the Heisenberg group G (H), with 6 
defined by 

(3.86) (q, P ) - e(V, P ')), q = (i+pTV, p = (i+pT'^-p) 

and with an inverse map (q',p') = C _1 (^(q,p)j given by q' = 2(1 + p)^ 1 q, p' = (1 +p) _1 (1 — p)- The 
Cayley transform is a conformal quaternionic contact diffeomorphism between the quaternionic Heisenberg 



(3.85) 



(-) 



[-dp + 2dq-q 



4 
4 
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group with its standard quaternionic contact structure and the sphere minus a point with its standard 
structure fj, [37]. In fact, by [37, Section 8.3] we have 

A-(e-^-A = ^e. 

where A = |1 +p | (1 + p)^ 1 is a unit quaternion and ij is defined in (3.84). 

Besides the Cayley transform, we shall need the generalization of the Euclidean inversion transformation 
to the qc setting. We recall that in [52] Koranyi introduced such an inversion and an analogue of the Kelvin 
transform on the Heisenberg group, which were later generalized in [22] and [21] to all groups of Heisenberg 
type. The inversion and Kelvin transform enjoy useful properties in the case of the four groups of Iwasawa 
type of which G (H) is a particular case. For our goals it is necessary to show that the inversion on G (H) 
is a qc-conformal map. In order to prove this fact we shall represent the inversion as the composition of two 
Cayley transforms, see [38, 44] where the seven dimensional case was used. Let P\ = (—1, 0) and P2 = (1, 0) 
be correspondingly the 'south' and 'north' poles of the unit sphere S 4n+3 = {\q\ 2 + \p\ 2 = 1} C H" x EL Let 
Ci and C2 be the corresponding Cayley transforms defined, respectively, on S 4n+3 \ {Pi} and S 4n+3 \ {i^}- 
Note that Ci was defined in (3.86), while C2 is given by {q,p) = 62 ( (<?',£>') 



(3.87) q = -(1-pT 1 q', p = (l-p'T 1 (I+P), (q', P ')£S An+3 \{P 2 } 

The inversion on the quaternionic Heisenberg group (with center the origin) with respect to the unit gauge 
sphere is the map 

(3.88) a = C 2 o e^ 1 : G (H) \ {(0, 0)} ->• G (H) \ {(0, 0)}. 



In particular, a = 62 Ci 1 is an involution on the group. A small calculation shows that a is given by the 
formula (in the Siegel model) 



* -1 

q = -p q, 



q* = -(M 2 - w ) 1 q, w * = -1 



or, equivalently, in the direct product model G (H) 

q* = 

Recalling (3.85) and (3.83) it follows 

a" 

(3.89) 



1 p 
9 = —rrr^Qp, p — — , (in the Siegel model) 
\P? \P\ 



<t* = -r-r, ; — ^rU&u, u = 77-, (in the product model) . 

M 4 + M 2 (M 4 + M 2 ) 1/2 

Thus, we have the following fundamental fact. 

Lemma 3.13. The inversion transformation (3.88) is a qc-conformal transformation on the quaternionic 
Heisenberg group. 

As usual, using the dilations and translations on the group, it is a simple matter to define an inversion 
with respect to any gauge ball. We recall that the gauge norm of a point (q,u>) £ G (H) is N(q,uj) = 
(|g| 4 + |w| 2 ) 1//4 which allows to define a distance on the group using the translation structure. 

After these preliminary facts, our next goal is to prove a version of Liouville's theorem in the case of the 
quaternionic Heisenberg group and the 3-Sasakian sphere equipped with their standard qc structures. In the 
Euclidean case Liouvillc [63] , [64] showed that every sufficiently smooth conformal map ( C 4 in fact) between 
two connected open sets of the Euclidean space R 3 is necessarily a restriction of a Mobius transformation. 
The latter is the group generated by translations, dilations and inversions of the extended Euclidean space 
obtained by adding an ideal point at infinity. Liouville's result generalizes easily to any dimension n > 3. 
Subsequently, Hartman [35] gave a proof requiring only C 1 smoothness of the conformal map, see also [67], 
[7], [49], [48] and [29] for other proofs. A CR version of Liouville's result can be found in [77] and [1]. Thus, 
a smooth CR diffcomorphism between two connected open subsets of the 2n + 1 dimensional sphere is the 
restriction of an element from the isometry group SU(n + 1, 1) of the unit ball equipped with the complex 
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hyperbolic metric. The proof of Alexander [1] relies on the extension property of a smooth CR map to 
a biholomorphism. Tanaka, see also [73], [13, 14] and [19], in his study of pscudo-conformal equivalence 
between analytic real hypersurfaces of complex space showed a more general result [77, Theorem 6] showing 
that any pseudo-conformal homeomorphism between connected open sets of the quadric 

r n 
i—l i=r+l 

is the restriction of a projective transformation of P ra (C). 

A new stage began with the introduction of quasiconformal maps, which imposed metric conditions on the 
maps, [32] and[75] and with the works of Mostow [66] and Pansu [70] . In particular, in [70] it was shown that 
every global 1-quasiconformal map on the sphere at infinity of each of the hyperbolic metrics is an isometry 
of the corresponding hyperbolic space. The local version of the Liouvillc's property for 1-quasiconformal 
map of class Cr on the Heisenberg group was settled in [53] by a reduction to the CR result. The optimal 
regularity question for quasiconformal maps in much greater generality, including the cases of all Iwasawa 
type groups, was proven later by Capogna [10], see also [78] and [11]. 

A closely related property is the so called rigidity property of quasiconformal or multicontact maps, also 
referred to as Liouvillc's property, but where the question is the finite dimensionality of the group of (locally 
defined) quasiconformal or multicontact maps, see [80], [74], [23], [68], [69], [65], [24], [56]. 

In Theorem 1.4 we give a local Liouville type property in the setting of a sufficiently smooth qc-conformal 
maps relying only on the qc geometry. 



3.9. Proof of Theorem 1.4. Since in the case £ = S An+3 there is nothing to prove we shall assume that 
S / S 4n+3 . We shall transfer the analysis to the quatcrnionic Heisenberg group using the Cayley transform, 
thereby reducing to the case of a qc-conformal transformation F : £ — ► G (H) between two domains of the 
quaternionic Heisenberg group such that 9 = F*Q — ^9 for some positive smooth function fi defined on 
the open set £. By its definition is a qc-Einstein structure, hence the proof of [37, Theorem 1.1] shows 
that fi satisfies a system of partial differential equations whose solution is a family of polynomial of degree 
four. In fact, with a small change of the parameters in [37, Theorem 1.1], \i has the form 



(3.90) ufau) = c 



g | 2 ) + \uj + uj + 2Img g| 



for some fixed (q a , w ) € G (H) and constants Co > and a G R. A small calculation using (3.90) and the 
Yamabe equation [37, (5.8)] shows Seals = 128n(n + 2)coct. Since 9 is qc-conformal to 9 via the map F, 
it follows that a = 0. Hence, F is a composition of a translation, cf. (3.82), followed by an inversion and a 
homothety, cf. Lemma 3.13. 

The above analysis implies that F is the restriction of an element of PSp(n + 1, 1). This completes the 
proof of Theorem 1.4. 

Similar to the Riemannian and CR cases Theorem 1.4 and a standard monodromy type argument, sec 
[54] and [76, Theorem VI. 1.6], show the validity of the next 

Theorem 3.14. If (M,ij) is a simply connected qc-conformally flat manifold of dimension An + 3, then 
there is a qc-conformal immersion $ : M — > S 4n+3 , where S 4n+3 is the 3-Sasakian unit sphere in the 
(n + 1) -dimensional quaternion space. 

3.10. Proof of Theorem 1.3. From Lemma 3.10, Lemma 3.11 and the classical Obata theorem it follows 
that (M, h) is isometric to the unit sphere in Euclidean space, i.e., there is a diffeomorphism i : M — > S 4n+3 
such that h = i*dx 2 , where dx 2 denotes the round metric on S in+3 which we take to be of constant 
Riemannian scalar curvature Scal h = (An + 3)(4n + 2). Thus, the curvature tensor R h of the Levi-Civita 
connection V' 1 of h is given by 



(3.91) 



R h {A, B, C, D) = h(B, C)h(A, D) - h(B, D)h(A, C). 
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The relation between the curvature tensors of the Levi-Civita and the Biquard connection [37, Corollary 4.13] 
or [44, Theorem 4.4.3] together with (3.91) yield 

(3.92) R(X, Y, Z, V) = g(Y, Z)g(X, V) - g(Y, V)g(X, Z) 



6 

Y [w.(y, Z)w a (X, V) - w,(X, Z)uj s (Y, V) - 2uj s (X, Y)lu s {Z, V) 



In the case T° = U — 0, S = 2, the formula for the qc-conformal curvature tensor given in [42, Proposition 4.2] 
reads 

1 r 3 i 

(3.93) W« C (X, Y, Z,V) = - R(X, Y, Z,V) + Y RQ.X, I.Y, Z, V) + g{X, Z)g(Y, V) - g(Y, Z)g(X, V) 



8=1 



Y ["s (X, Z)uj s (Y, V) - w s (Y, Z)uj s (X, V) 



s=l 

With a small calculation we see from (3.93), taking into account (3.92), that the qc-conformal curvature 
tensor vanishes, W qc = and (M, g, r), Q) is locally qc-conformal to the sphere due to [42, Theorem 1.3]. 
Hence, taking into account Theorem 3.14, (M,g, r], Q) is qc-conformal to S 4n+3 , i.e., ?y = K^F*fj for some 
diffeomorphsm F : M —> S 4n+3 . Comparing the metrics we obtain n = 1 which shows that M is qc- 
homothctic to the 3-Sasakian unit sphere in the (n + l)-dimcnsional quaternion space. This completes the 
proof of Theorem 1.3. The proof of Theorem 1.2 follows as already noted after the statement of Theorem 1.3. 

4. Appendix 

4.1. The P— form. Let (M, g, Q) be a compact quaternionic contact manifold of dimension An + 3 and / a 
smooth function on M. We recall the notion of a P-function introduced in [41] 

Definition 4.1. /[41]/ 

a) For a fixed f we define a one form P = Pj = P[f] on M , which we call the P—form of f , by the following 
equation 



P f (X) = V 3 f(X,e b ,e b ) + Yv 3 f(ItX,e b ,I t e b ) - 4nSdf(X) + 4nT°(X,Vf) 



in(n-2) 
n-l 



U(X,Wf). 



b) The P— function of f is the function Pf(Vf). 

c) The C— operator is the fourth- order differential operator on M (independent of f!) defined by 

Cf=-V*Pf = (V ea Pf){e a ). 

d) We say that the P— function of f is non-negative if its integral exists and is non-positive 



(4.1) 



/ f-CfVol v = - [ P f (Vf)Vol v >0 
Jm Jm 



If (4.1) holds for any smooth function of compact support we say that the C '— operator is non-negative. 
The 5p(n)5p(l)-invariant decomposition of the horizontal Hessian V 2 / are given by 

1 



(v 2 /) [3] (x,y) = - V 2 f{X,Y) + Y v 2 f{i s x,i s Y) 



(4.2) 



s=l 
3 



(V 2 /) M] (X, Y) = - 3V 2 /(X, Y)~Y V 2 f(LX, I S Y) 



Let (V 2 /)[3][oj be the trace-free part of the 3-component of the horizontal Hessian, 



(4.3) 



(S7 2 f) [m (X,Y) = (V 2 /) [3] (X,r) + —Afg(X,Y). 
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The next local formula, established in [41], 

(4-4) (V eo (V 2 /)[3][0])(ea,X) 



n- 1 
An 



implies the non-negativity of the C— operator on a compact qc manifold of dimension at least eleven [41, 
Theorem 3.3]. Indeed, using (4.4) we have 



(4.5) 



n-1 
An 



f.CfVoL = - 



AI 



n-1 

An 



A I 



P f (Vf)Vol v 



A I 



|(V7)[3][0]| 2 ^ r „ 



after using an integration by parts and the orthogonality of the components of the horizontal Hessian. 

4.2. A new proof of Theorem 1.1. Here we use the non-negativity of the P-function P(V/) of a smooth 
function / established in [41, Theorem 3.3] to give a new proof of Theorem 1.1. 

Proof. Let / be an eigenfunction of the sub-Laplacian with eigenvalue A, i.e., we assume that (2.14) holds. 
An integration by parts gives 



(4.6) 



/ (A/) 2 Vol v = A / fAf Voir, = A / |Vjf Voir, 

JM JM JM 



We recall the qc-Bochncr identity [40, (4.1)]. Applying the first equality in (2.10), (2.9) and the properties 
of the torsion, (2.8), we can write the qc-Bochner formula [40, (4.1)] in the form 



(4.7) ^ A|V/I = |V 71 " 9 ( y ( A ^)' V /) + 2 ( n + 2 )SI V /r + 2(n + 2)T°(V/, V/) 



+2(2n + 2)U(Vf, V/) + 4 £ V 2 /(6, W). 



s=l 



One of the key identities which relates the P-function and the qc-Bochner formula (4.7) is given by the 
following equation [41, Lemma 3.2] 



(4.8) 



' Vv 2 /(6,/ s v/)H 



A I 



C/(V/,V/) 



Vol r . 



An integration of (4.7) over the compact M, followed by a substitution of (2.14) and (4.8) in the obtained 
integral equality, and then a use of the divergence formula (2.16) give 



(4.9) 



A I 



|V 2 /| 2 - A|V/| 2 + 2n5|V/| 2 + 2(n + 2)T°(V/, V/) + ^ \ 1} U(Vf, V/) 



n- 1 
n n 



V0l v . 



The latter formula can be written in the form 
= 

2n + 1 



(4.10) 



2(n + 2) 



f { | V 2 /| 2 - A| V/| 2 - 5| V/| 2 + T°(V/, V/) - UpfV/, V/) 

2( „ + 2)S ,V/P + 4 " 2 + "" +12 T°(V/, V/, + t^fV'V /. V/)] 
2n + l (n — l)(2n+l) J 

--p„(v/)--(A/) 2 } y z„. 

Now we invoke the next integral identity proved in [40, Lemma 3.4] 

3 „ 3 3 



(4.11) [ J2\7 2 f({ s ,I s Vf)Vol v = - [ k^(d/(£ s )) 2 + ^T(£ s ,/ s V/,V/) 

J M s=1 J M 1 S =l s=l 



Vol v . 
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From (4.11) and (4.8) it follows the equality 

(4.12) / - s\ v/| 2 + r°(v/, v/) - 2{n ~? U(Vf, V/) 



M 



n — 1 



Vol n = 



i-^[ 5 (V 2 /^.s)] 2 }^ 



A substitution of (4.12) in (4.10) yields 

1 



(4.13) 0=/ {|V 2 /| 2 --[(A/) 2 + ^[.9(V 2 /,^)] 



I M 

2n + 1 r 



2S\Vf\ 



An 2 + 14n+ 12, 
(2n+ l)(n + 2)' 



4(n + 2)(2n - 1) 



r°(v/, v/) + in+ "; M v/,v/) - ^|v/| 2 ] } ^ 

(n — l)(2n + 1) n i > 



The equality (4.13), the Lichnerowicz type assumption (3.4) and (4.6) imply the inequality 



(4.14) 0>/ {|V 2 / | 2 -^[(A/) 2 + ^[ 5 (V 2 / , Ws )] 2 ]-£p4V/) + ^( 



2n + 1 / k A 



n + 2 n 



|V/| 2 } Voir,. 



Note that in the proof of (4.14) we supposed implicitly that n > 1. But it works also for n = 1, when [7 = 
trivially, we have only to remove the torsion tensor {/ (cf. [41]). 

Using that | 2^777^} ^ s an orthonormal set in ^[-i] we have 



(4.15) 



l(v 2 /) [ - 1] l 2 >^E[.g(v 2 /,^)] ; 



while a projection on {^T^ff} gi yes 



(4.16) 



|(V 2 /) [3] | 2 >-(A/) 2 



Next, using the Sp(n)Sp(l)-invariant orthogonal decomposition (4.2) of horizontal Hessian and the estimates 
(4.15) and (4.16), we obtain the inequality 



(4.17) 



V 2 /| 2 >i-[(A/) 2 + E[.g(V 2 /, Ws ^ 



in 



s=l 



Finally, using the non- negativity of the P- function for n > 1, see (4.5), we obtain from (4.14) the desired 
estimate 

n , 



A > 



n + 2 



□ 



Corollary 4.2. If the case of equality in Theorem 1.2 holds, i.e., we have 



A = 



n + 2 



A/ = — -kof, 
n + 2 



then the horizontal Hessian of the eigenfunction f is given by (1.4). 

Proof. The result follows from (4.13), (3.4) and (4.17) which asserts that the equalities in (4.15) and (4.16) 
must hold, which imply (1.4) . □ 
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